MAY 1<± MI KB No. 3J30 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



WARTIME REPORT 



Restricted Bulletin 3J30 



USE OF STAGNATION TEMPERATURE IN CALCULATING RATE OF 
HEAT TRANSFER IN AIRCRAFT HEAT EXCHANGERS 
By George P. Wood 



Langley Memorial Aeronautical Laboratory 
Langley Field, Va. 



MAC A 



!! ' WASHINGTON 

safe 

li 

[1 NACA WARTIME REPORTS are reprints of papers originally issued to provide rapid distribution of 
; advance research results to an authorized group requiring them for the war effort. They were pre- 
viously held under a security status but are now unclassified. Some of these reports were not tech- 
nically edited. All have been reproduced without change in order to expedite general distribution. 



L " LANGLEY MEMORIAL AliRONAUTICAl 

LABORATORY 
Langley Field, V«l 

li. — ^ - 



1176 01354 2643 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



RESTRICTED BULLETIN 



USE OP STAGNATION TEMPERATURE IN CALCULATING RATE OP 
HEAT TRANSFER IN AIRCRAFT HEAT EXCHANGERS 
By George P. Wood 



SUMMARY 



Theoretical and experimental investigations of the 
effect of friotional heat on the rate of heat transmis- 
sion at high fluid velocities are briefly reviewed. On 
the basis of these investigations, calculations are made 
which show that the use of the stagnation temperature of 
the cooling air as the effective temperature for heat 
transfer in an aircraft heat exchanger is sufficiently 
accurate. 



INTRODUCTION 



No new analyses or experiments are given in the 
present paper. Its purpose is to point out that, on the 
basis of existing information, the proper cooling-air 
temperature for use in aircraf t-heat-exchanger calcula- 
tions is, to b sufficient degree of accuracy, the stag- 
nation temperature of the cooling air. 

At the present time some uncertainty exists con- 
cerning the proper cooling-air temperature to be used 
in calculations of rate of heat transfer when the 
velocity of the cooling air across the heat-transfer 
surface Is high, or when the velocity is low but the 
air has been slowed down from a high-velocity air stream. 
Conventional heat-transf rr coefficients generally have 
been determined from measurements made at rather low air 
velocities. At high velocities, the heat generated In 
the boundary layer by friction may bo an important fac- 
tor in the rate at which heat is transferred across the 
boundary layer from the surface to the fluid. The effect 
of the frictional heat, which is negligible at low air 
velocities, is not included in the -conventional low- 
velocity coefficients of heat transfer or in the 
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conventional definition of the temperature difference 
between the surface and the fluid. The problem of 
finding the equation for rate of heat transfer that pro- 
perly takes into account the conditions in the boundary 
layer for high fluid velocities has been undertaken by a 
number of investigators. Because the mathematics of the 
problem is complex and simplifying assumptions must be 
made in order to obtain solution's, rigorous results for 
all types of flow have not yet been obtained. Because 
many different equations have been derived, some confu- 
sion exists. Actually, when typical numerical values 
are substituted in the various equations, the results, 
though differing slightly in numerical value, are all of 
the same order of magnitude. The present paper shows 
that, for aircraft heat exchangers, a rigorous equation, 
for high-velocity heat transfer is not necessary but 
that, on the basis of the investigations reviewed, suf- 
ficient accuracy results from the use of the stagnation 
temperature as the effective temperature of the cooling 
air. 



SYMBOLS 



specific heat of fluid at constant pressure, 
Btu per pound per °P 



f 



friction factor, dimenslonless 
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acceleration due to gravity, feet per second 
per second 



h 



coefficient of heat transfer for low fluid 
velocities, Btu per second per square foot 
per °P 



H 



J 



rate of heat transfer, Btu per second 
mechanical equivalent of heat (778 ft-lb/Btu) 



k 



thermal conductivity of fluid, Btu per second 
per square foot per °F per foot 



P 



rate of doing v/ork against friction, foot-pounds 
per second 



Pr Prandtl number, dimensionless, (o p ng/k) 

S surface area for heat transfer, square feet 

T temperature,, °F, . 

AT temperature difference, °F 

T Q free-stream temperature of fluid, °P 

(AT.) ad adiabatic temperature change due to compression, 

Op 

T eff fluid temperature effective for heat transfer, °P 

(AT) eff difference between T eff and T f , 

If actual temperature of fluid outside boundary 
layer, °P 

T__ temperature of fluid at heat-exchanger entrance, 
en 0 p 

T ex temperature of fluid at heat-exchanger exit, °P 

T a stagnation temperature of fluid, °P 

T a stagnation temperature of fluid at heat- exchanger 
en entrance, °P 

T. stagnation temperature of fluid at heat-exchanger 
ex exit, °F 

T w temperature of wall, °P 

V velocity of fluid, feet per second 

V Q free-stream velocity of fluid, feet per second 

V velocity of fluid at heat-exchanger entrance, 

feet per second 

V ex velocity of fluid at heat-exchanger exit, feet 
per second 

W weight rate of flow of fluid, pounds per second 

y normal distanoe from flat plate, feet 
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coefficient of viscosity of fluid, slugs per 
foot-second 



P 



mass density of fluid, slugs per cubic foot 



T 



frictional resistance per unit surface area, 
pounds per square foot 



DISCUSSION 



Consider a flat plate with its axis along the flow 
direction. The rate at which heat is transferred between 
the plate and the fluid Is determined by the temperature 
gradient at the plate in the boundary layer of the fluid 
and by the thermal conductivity of the fluid; that is, 



The method of determining, from the differential equa- 
tions of the boundary layer, the effect of the frictional 
heat on the temperature gradient need not be discussed 
in the present paper. An Idea, however, of the reason 
that the frictional heat may be an important factor at 
high speeds can be had from the following simple analysis: 

The frictional resistance per unit surface area Is 
the product of friction factor and dynamic pressure 



The rate at which frictional work is done per unit sur- 
face area is 
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The rate of heat transfer per unit surface area ±s. 



H 
3 



h AT 



The ratio of work done to heat transferred Is then 



P 
H 



tPV* 
2h AT 



By using the relation discussed In reference 1 for cor- 
relating friction and heat transfer 



the ratio between frictlonal work and transferred heat 
is obtained as 



The ratio of the heat generated In overcoming friction 
to the heat transferred between surface and fluid is 
thus proportional to the square of the fluid velocity. 
It Is therefore to be expected that the effect of the 
frictlonal heat Is most pronounced at high velocities. 

The problem of obtaining a rigorous expression for 
rate of heat transfer with high fluid velocities has 
been the subject of a number of Investigations. Some of 
the results that have been obtained are briefly discussed. 

For laminar flow, analyses (references 2 and 3) show 
that the proper temperature difference for use in calcu- 
lating rate of heat transfer for high velocities and 
with use of the usual low-velocity coefficients Is the 
difference between the wall temperature T w and the ef- 
fective - temperature T e ff, where T e ff Is a function of 



2h = 



CpgpVf 




P _ V 2 (Pr ) g / 3 
H " c p g AT 



(1) 
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the Prandtl- numben of the fluid, the fluid temperature 
outside the boundary layer Tf , and the increase In 
temperature that would result from complete adiabatic 
arrest of the flow outside the boundary layer (AT) ad . 
The quantity 



( * T >eff = T eff- T f 



can be found from figure 1, which- is taken from figure 
2 of reference 2. 



(AT) 



eff 



(AT) 



ad 




1 4 10 40 100 400 10C0 
Prandtl number, Pr 

Figure 1.- Ratio (AT) eJ?f /( AT) ad as a function of Prandtl 
nunber for laminar flow. (Prom fig. 2 of reference 2.) 



For air, the Prandtl number of which. is about 0.75, 
the value of (AT) eff /(AT) ad is approximately 0.85. The 

rate of heat transfer between an element of wall area S 
and air in. laminar flow is then caloulat«jd by the equa- 
tion 



H 

S 



= h( k T w -T ef f ) 

= h[T w -T f -(AT) eff ] 
= h[T w -T f -0.85(AT) ad J 
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For a given fluid, (AT) d (In ®P) Is a function- of only 
velocity. For air, 

(AT)- ad = (0.832xlO- 4 )V 2 



That the temperature T 6 ff produced at the wall by 
friction should be a function of the Prandtl number Pr 
of the fluid seems reasonable because the rate at which 
heat is generated In the boundary layer is proportional 
to the coefficient of viscosity of the fluid n» Inasmuch 
as 



T =3 



and because the rate at which heat is conducted ^.cross 
the boundary layer is proportional to the thermal con-, 
ductivity of the fluid k. 

Experiments made with air flowing past unheated 
surfaces, which are reported in references 4 and 5-, have 
corroborated the relation 



(AT) eff = 0,85(AT) ad 



for laminar flow. (See fig. 2 which Is taken from 
fig. 4 of reference 4.) 



(AT) 



eff 



(AT) 



ad 



Laminar 
.90 



I 



Turbulent 




.80 



10 5 2 4 6 10 6 2 4 
Reynolds number 



Figure 2.- Ratio (AT) ff /(AT) ad .as a function of 
Reynolds number. Axial wire In flow of air through 
nozzle. (From fig. 4 of reference 4.) 
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For turbulent flow, uncertainty ^concerning the exact 
distribution of. velocity In the laminar sublayer at the 
fluid boundary- and concerning conditions at the Interface 
between. the laminar sublayer and the buffer layer, to- 
gether with the difficulty of 'obtaining a solution to the 
differential equations of the boundary layer, has led to 
somewhat different results In different Investigations. 
Reference 6 gives as the effective fluid temperature T e ff 
the full stagnation temperature, with rate of heat 
transfer calculated by the equation 

| =\ h (V T eff) 

= .h[T w -T- f -(AT) ad ] 

The result, however, is not rigorous because the laminar 
sublayer is not taken into account. In reference 7 
another relation is derived but, according to a statement 
in reference 8, the analysis Is based on relatively rough 
assumptions. In reference 8 the analysis takes into ac- 
count the laminar sublayei , but both the original' and the 
revised analyses given in reference 8. appear to be bas>ed 
on assumptions that may not be wholly correct.- The re- 
vised analysis also appears to contain algebraic mistakes. 
In reference 9 the problem is treated only for fluids 
that have a Prandtl number of unity. For such fluids, 
T eff ls founcl to De the full stagnation temperature and 




Many other references might be mentioned - for example, 
about 20 papers related to the subject are referred to 
in reference 10. None of them, however, seems to be 
entirely unquestionable, but all derive answers that, for 
air, are of the same order of magnitude. As will be 
shown later, however, the rigorous solution of the pro- 
blem is not necessary for practicable beat-exchanger cal- 
culations, and only the order of magnitude of the answer 
is needed. 
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The only experimental determination of T 9 ff for 
high-velocity turbulent flow appears to be some measure- 
ments made with air and reported in reference 4. These 
- measurements consisted of the. determination of the tem- 
perature assumed by a wire placed axial ly in a. nozzle. 
The results are shown in figure 2. Figure 2 shows that 
the ratio (AT) eff /(AT) ad is slightly higher for turbu- 
lent flow than for laminar flow and # at least for the 
range of Reynolds number shown, is of the order of 0.88. 

For the entrance region of a heat exchanger, in 
which the velocity profile in the fluid is strongly de- 
pendent upon the distance from the leading edge, no 
analysis has been made. 

Equation (1) shows that the pumping power for a 
heat-transfer surface is proportional to the square of 
the fluid Velocity. That the fluid velocity should be 
kept low in heat exchangers is generally realized. The 
discussion in the present paper of the correct fluid 
temperature for heat-transfer calculations with high 
fluid velocities is therefore primarily of only academic 
interest. (An exception, however, is the relatively 
rare case of heat transfer from an airplane wing.) In 
most heat exchangers, the air velocity is low enough 
that the effect of frictlonal generation of heat on the 
air temperature is unimportant. On the other hand, the 
airplane velocity is generally high enough that the 
adiabatlc rise in temperature of the air between the 
free stream and the entrance to the exchanger must be 
taken into account. Inasmuch as the indications are, 
as has been previously discussed, that the effective 
temperature rise is given by practically the same 
equation as the compressional stagnation temperature 
rise and because the actual compressional rise generally 
is almost complete stagnation rise, calculations of rate 
of heat transfer in a heat exchanger will usually be 
sufficiently accurate if the temperature of the cooling 
air at the entrance to the heat exchanger is taken to be 
the full stagnation temperature of the free-stream air 



Ten = T a = T Q +(AT) ad 
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TABLE I 

COMPARISON OP ACTUAL, STAGNATION, AND EFFECTIVE AIR 
TEMPERATURES FOR VARIOUS VALUES OF V Qn AND FOR 

V Q ■= 500 mph 



v en 


T-, - T 
f ~ o 


T - T 

A s x o « 


T eff " T o 


(mph) 


(°F) 


<°F) 


<°F) 


400 


16.2 


45 


42.1 
■ 


300 


28.8 


45 


43.4 


200 


37.8 


45 


44.3 


100 


43.2 


46 


44.8 



Table I permits a comparison of the Values of the 
actual Inlet temperature Tf and the stagnation tempera- 
ture T fl with the approximate values of tne effective 
temperature T ef ..p. for various entrance velocities and 

for a free- stream velocity of 500 miles per hour. The 
values shown in the table were calculated as follows: 

For a velocity of 200 miles per ho-ur in the. heat- 
exchanger .entrance,- the actual fluid temperature, in. the 
exchanger- entrance '.in °F- is 

T = T+(l. 8x10-4) [(500) 8 - (200) 2 1 
en o " ' ■_ _j 

- T 0 +37.8 

Stagnation temperature in °F ia 

T s = T o +{l.8xl0- 4 ) (500) 2 
= T. 0 +45 
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The temperature that is effective for heat transfer 
T eff ln °& ls » according to the data of figure 2, of 
the order of 

T eff = T Q + (l.8xl0-*) [( 500) 2 -( 200) 2 +0.9 (200) 2| 
= Tq+44.3 



A comparison of the stagnation temperatures with 
the effective temperatures shown in the last two columns 
of table I shows that, at all cooling-air velocities, 
the difference between the two temperatures is insigni- 
ficant. It therefore appears that use of full stagna- 
tion temperature as the temperature effective for 'heat 
transfer is sufficiently accurate for the calculation of 
heat transfer in aircraft heat exchanger?. 

Two further interesting facts can be mentioned. In 
calculating heat transfer at high altitudes, it is more 
desirable to use the exact heat-balance equation 

H = W°p(Tex-Ten)-^(V e x 2 - v en 2 ) < 2 > 



than the inexact but more generally used .equation 
H = Wc p( T ex- T en) 



At low altitudes, the last term in equation (2) is usual- 
ly negligible but, at high altitudes, may be important. 
If the effective temperature is taken as the stagnation 
temperature 



T er _ T eff ~ T s 



calculations in which equation (2) is used are greatly 
simplified. The simplification results from the fact 
that equation (2) can be written 
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H 



P(^ s ex' Ts 



en 



It should also be pointed out that use of the stag- 
nation temperatures of the cooling air and the engine air 
instead of the actual fluid temperatures in calculations 
for cross-flow intercoolers in no way invalidates 
Nusselt's determination of the mean temperature difference 
in cross flow (reference 11). Nusselt's calculation of 
mean temperature difference is based on the actual tem- 
peratures of the fluids and on the approximate heat- 
balance equation (3) . An inspection, however, of his 
analysis shows that if, in the analysis, actual tempera- 
tures are replaced by stagnation temperatures and the 
approximate heat-balance equation (3) is replaced by 
the exact equation (2), the analysis proceeds exactly 
as riven by Nusaelt and the same numerical results are 
obtained. 



A brief review has been given of the investigations 
of the effect of boundary- layer heat on rate of heat 
transmission. Calculations based on the results of 
these investigations showed that sufficient accuracy can 
be obtained in calculations of the rate of heat transfer 
in aircraft heat exchangers by using the stagnation tem- 
perature of the cooling air as its effective temperature. 



CONCLUDING RTSMARKS 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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